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This paper reviews the theoretical approaches and existing models of the solar wind 
interaction with the Local Interstellar Cloud (LIC). Models discussed take into account 
the multi-component nature of the solar wind and local interstellar medium. Basic results 
of the modeling and their possible applications to interpretation of space experiments are 
summarized. Open questions of global modeling of the solar wind/LIC interaction and 
future perspectives are discussed. 

1. INTRODUCTION 

The Local Interstellar Cloud (LIC) is a cloud of partly ionized plasma surrounding 
the Solar System. The plasma component of LIC interacts with the solar wind plasma 
and forms the heliospheric interface (Figure 1). The heliospheric interface is a complex 
structure, where the solar wind and interstellar plasma, interplanetary and interstellar 
magnetic fields, interstellar atoms of hydrogen, galactic and anomalous cosmic rays (GCRs 
and ACRs) and pickup ions play prominent roles. 

Although a space mission into the Local Interstellar Cloud is becoming now more real- 
isable, there are no yet direct observations inside the heliospheric interface. Therefore, at 
the present time the heliospheric interface structure and local interstellar parameters can 
be derived only from remote experiments and measurements. Currently, backscattered 
solar Ly-a radiation, pickup ions, anomalous cosmic rays, and Voyager measurements of 
distant solar wind are the major sources of information on the heliospheric interface struc- 
ture and position of the termination shock ffl. kHz emission detected by Voyager can put 
some constraints. Recently, it was shown that study of Ly-a absorptions toward nearby 
stars can serve as remote diagnostics of the heliospheric interface and, in particular, the 
hydrogen wall around the heliopause (e.g., @ - H). In the foreseeable future, remote 
diagnostics will be also possible with images of heliospheric energetic neutrals (ENAs) . 
To reconstruct the structure of the interface and physical processes inside the interface 
on the basis of remote observations, a theoretical model should be employed. 

Theoretical studies of the heliospheric interface were performed over more than four 
decades after pioneering papers by Parker |7[] and Baranov et al. ||. However, a com- 
plete theoretical model of the heliospheric interface has not been constructed yet. The 
difficulty in doing this is connected with the multi-component nature of both the LIC 
and the solar wind. The LIC consists of at least five components: plasma (electrons and 
protons), hydrogen atoms, interstellar magnetic field, galactic cosmic rays, and interstel- 
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Figure 1. The heliospheric interface is the region of the solar wind interaction with 
LIC. The heliopause is a contact discontinuity, which separates the plasma wind from 
interstellar plasmas. The termination shock decelerates the supersonic solar wind. The 
bow shock may also exist in the interstellar medium. The heliospheric interface can be 
divided into four regions with significantly different plasma properties: 1) supersonic solar 
wind; 2) subsonic solar wind in the region between the heliopause and termination shock; 
3) disturbed interstellar plasma region (or "pile- up" region) around the heliopause; 4) 
undisturbed interstellar medium. 



lar dust. The heliospheric plasma consists of original solar particles (protons, electrons, 
alpha particles, etc.), pickup ions and the anomalous cosmic ray component. The pickup 
ion component is a result of ionization of those interstellar H atoms that penetrate into 
the heliosphere through the heliospheric interface. A part of the pickup ions is accelerated 
to high energies of ACRs. ACRs may also modify the plasma flow upstream of the ter- 
mination shock and in the heliosheath. Spectra of ACRs can serve as remote diagnostics 
of the termination shock. For a recent review on ACRs see ||. 

To construct a theoretical model of the heliospheric interface, one needs to choose a spe- 
cific approach for each interstellar and solar wind component. Interstellar and solar wind 
protons and electrons can probably be described as fluids. At the same time interstellar 
H atom flow requires kinetic description. For pickup ion and cosmic ray components, 
the kinetic approach is also required. However, for interpretations that are not directly 
connected to pickup ions and ACRs, a cruder model can be used. 
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Table 1 

Number Densities and Pressures of Solar Wind Components 



Component 



4-5 AU 
Number Density 
cm -3 



Pressure 
eV/cm~ 3 



80 AU 
Number Density 
cm -3 



Pressure 
eV/cm~ 3 



Original solar 
wind protons 
Pickup ions 
Anomalous 
cosmic rays 



0.2-0.4 
5.1 • 10~ 4 



2.-4. (thermal) 
^ 200 (dynamic) 
0.5 



(7-14) ■ 10 



2 • 10" 4 



-4 



10-3 - 10" 4 
0.5 — 1. (dynamic) 
~ 0.15 

0.01 - 0.1 



This paper focuses on the models of the global heliospheric interface structure. Under 
global models I understand those models that study the whole interaction region, including 
the termination shock, the heliopause and possible bow shock. In this sense, this paper 
should not be considered as a complete review of progress in the field. Many different 
approaches were used to look into different aspects of the solar wind interaction with 
LIC connecting with pickup ion transport and acceleration, with the termination shock 
structure under influence of ACRs and pickup ions. For more complete overview see recent 
reviews |T|,@. 

The structure of the paper is the following: The next section briefly describes our 
current knowledge of the local interstellar and solar wind parameters. Section 3 discusses 
theoretical approaches to be used for the interstellar and solar wind components. Section 
4 gives an overview of heliospheric interface models. Section 5 describes basic results 
of the Baranov-Malama model of the heliospheric interface and its future developments. 
In section 6 we demonstrate possible analyses of space experiments on the basis of a 
theoretical model of the heliospheric interface. Section 7 underlines current problems in 
the modeling of the global heliosphere and discusses future perspectives. 

2. BRIEF SUMMARY OF OBSERVATIONAL KNOWLEDGE 

Choice of an adequate theoretical model of the heliospheric interface depends on bound- 
ary conditions, i.e. on undisturbed solar wind and interstellar parameters. 

2.1. Solar wind observations 

At the Earth's orbit the flux of interstellar atoms is quite small, and the solar wind 
can be considered undisturbed. Measurements of pickup ions and ACRs also show that 
these components do not have dynamical influences on the original solar wind particles 
at the Earth's orbit. Therefore, solar wind parameters at the Earth's orbit can be taken 
as inner boundary conditions. 

It has been shown by many authors that pickup and ACR components dynamically 
influence the solar wind at large heliocentric distances. Observable evidence of such 
influence is, for example, deceleration of the solar wind detected by Voyager fi~I| . Table [I] 
presents estimates of dynamic importance of the heliospheric plasma components at small 
and large heliocentric distances. The table shows that pickup ion thermal pressure can 
be up to 30-50 % of the dynamic pressure of solar wind. 
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Table 2 

Local Interstellar Parameters 

Parameter 

Sun/LIC relative velocity 

Local interstellar temperature 

LIC H atoms number density 

LIC proton number density 

Local Interstellar magnetic field 

Pressure of low energetic part of cosmic rays 



Direct measurements/estimations 

25.3 ± 0.4 km s" 1 (direct He atoms x ) 

25.7 km s" 1 (Doppler-shifted 

absorption lines 2 ) 

7000 ± 600 K (direct He atoms x ) 

6700 K (absorption lines 2 ) 

0.2 ± 0.05 cm -3 (estimate based on 

pickup ion observations 3 ) 

0.03 - 0.1 cm -3 (estimate based on 

pickup ion observations 3 ) 

Magnitude: 2-4 fiG 

Direction: unknown 

~0.2 eV cm" 3 



2.2. Interstellar parameters 

Local interstellar temperature and velocity can be inferred from direct measurements of 
interstellar atoms of helium by Ulysses / GAS instrument . Atoms of interstellar helium 
penetrate the heliospheric interface undisturbed, because of the small strength of their 
coupling with interstellar and solar wind protons. Indeed, due to small cross sections of 
elastic collisions and charge exchange with protons, the mean free path of these atoms is 
larger than the heliospheric interface. Independently, the velocity and temperature in the 
Local Interstellar Cloud can be deduced from analysis of absorption features in the stellar 



spectra |12| . However, this method provides mean values along the line of sight in the 
LIC. A comparison of local interstellar temperatures and velocities derived from stellar 
absorption with those derived from direct measurements of interstellar helium shows quite 
good agreement (see Table 2). 

Other local parameters of the interstellar medium, such as interstellar H atom and elec- 
tron number densities, and strength and direction of the interstellar magnetic field, are not 
well known. In the models they can be considered as free parameters. However, measure- 
ments of interstellar H atoms and their derivatives as pickup ions and ACRs provide im- 
portant constraints on local interstellar densities and total pressure. The neutral H density 
in the inner heliosphere depends on filtration the neutral H atoms in the heliospheric inter- 
face due to charge exchange. Since interstellar He is not perturbed in the interface, local 
interstellar number density of H atoms can be estimated from the neutral hydrogen to the 
neutral helium ratio in the LIC, R(HI / Hel)nc'- niic(HI) = R(HI / HeI)Lic n Lic(HeI). 
The neutral He number density in the heliosphere has been recently determined to be very 
likely around 0.013 — 0.018 cm" 3 (JTB| - |T3|). Interstellar ratio HI/Hel is likely in the 
range of 10-14. Therefore, expected interstellar H atom number densities are in the range 
of 0.13 — 0.25 cm" 3 . It was shown by modeling [16] , [17]] that the filtration factor, which 



is the ratio of neutral H density inside and outside the heliosphere, is a function of in- 
terstellar plasma number density. Therefore, the number density of interstellar protons 
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(electrons) can be estimated from this filtration factor ||12|| . Independently, the electron 
number density in the LIC can be estimated from abundances ratios of ions of different 
ionization states 
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Note that there are other methods to estimate interstellar H atom density inside the 
heliosphere, based on their influence on the distant solar wind [11 or from ACR spectra 
E7| | . Recent estimates of the location of the heliospheric termination shock using transient 
decreases of cosmic rays observed by Voyager 1 and 2 also provide constraints on the 
local interstellar parameters [18|. However, simultaneous analysis of different types of 



observational constraints has not been done yet. Theoretical models should be employed 
to make such analysis. Table |2| presents a summary of our knowledge of local interstellar 
parameters. Using these parameters, we estimate local pressures of different interstellar 
components (Table 3). All pressures have the same order of magnitude. This means that 
theoretical models should not neglect any of these interstellar components. Dynamical 
pressure of interstellar H atoms is larger than all other pressures. A part of H atoms, ACRs 
and GCRs penetrate into the heliosphere, which makes their real dynamical influence on 
the heliospheric plasma interface difficult to estimate. 



Table 3 

Local Pressures of Interstellar Components 



Component 


Pressure estimation, dyn cm 2 


Interstellar plasma component 




Thermal pressure 


(0.6-2.0) • 10~ 13 


Dynamic pressure 


(1.5-6) • 10~ 13 


H atoms 




Thermal pressure 


(0.6-2.0)- 10~ 13 


Dynamic pressure 


(4.0 - 9.0) • 10- 13 


Interstellar magnetic field 


(1.0-5.0) ■ 10" 13 


Low energy part of GCR 


(1.0-5.0) ■ 10" 13 



3. OVERVIEW OF THEORETICAL APPROACHES 

In this section we consider theoretical approaches for components involved in the dy- 
namical processes in the heliospheric interface. 

Generally, any gas can be described on a kinetic or a hydrodynamic level. In the ki- 
netic approach, macroscopic parameters of a gas of s-particles (or, briefly, s-gas) can be 
expressed through integrals of velocity distribution function f s (f,w,t): n s = J f s dw, V s = 
(Jwf s dw)/n s , P sAj = m s S{w i -V 8t i)(w j -_V s>j )f s dw, q s = 0.5m s J(w - V s ) 2 (w - V s )f s dw, 
where n s is the number density of s-gas, V s is the bulk velocity of s-gas, P s>i j are compo- 
nents of the stress tensor P s , q s is the thermal flux vector, m s is the mass of individual 
s-particle. In the hydrodynamic approach, some assumptions should be made to specify 
the stress tensor P s , and the thermal flux vector, q s to make hydrodynamic system closed. 
For example, these values can be calculated by the Chapman-Enskog method, assuming 



6 



Kn = l/L << 1, where I and L are the mean free path of the particles and character- 
istic size of the problem, respectively. The zero approximation of the Chapman-Enskog 
method gives local Maxwellian distribution, and the gas can be considered as an ideal 
gas, where the stress tensor reduces to scalar pressure P and q — 0. 

3.1. H atoms 

Interstellar atoms of hydrogen form the most abundant component in the circumsolar 
local interstellar medium (see, Table 2). These atoms penetrate deep into the heliosphere 
and interact with interstellar and solar wind plasma protons. The cross sections of elastic 
H-H, H-p collisions are negligible as compared with the charge exchange cross section 
Charge exchange with solar wind/interstellar protons determines the properties of 
the H atom gas in the interface. Atoms, newly created by charge exchange, have the 
local properties of protons. Since plasma properties are different in the four regions 
of the heliospheric interface shown in Figure 1, the H atoms can be separated into four 
populations, each having significantly different properties. The strength of H atom-proton 
coupling can be estimated through the calculation of mean free path of H atoms in plasma. 
Generally, the mean free path (with respect to the momentum transfer) of s-particle in 
t-gas can be calculated by the formula: / = m s w^/(5M st /5t). Here, w s is the individual 
velocity of s-particle, and 6M st /5t is individual s-particle momentum transfer rate in t-gas. 

Table ^ shows the mean free paths of H atoms with respect to charge exchange with 
protons. The mean free paths are calculated for typical atoms of different populations at 
different regions of the interface in the upwind direction. For every population of H atoms, 
there is at least one region in the interface where the Knudsen number Kn ~ 0.5 — 1.0. 
Therefore, the kinetic Boltzmann approach must be used to describe interstellar atoms in 
the heliospheric interface. 



Table 4 

Mean free paths of H-atoms in the heliospheric interface with respect to charge exchange 
with protons, in AU. 



Population 


At TS 


At HP 


Between HP and BS 


LISM 


4 (primary interstellar) 


150 


100 


110 


870 


3 (secondary interstellar) 


66 


40 


58 


190 


2 (atoms originating in the heliosheath) 


830 


200 


110 


200 


1 (neutralized solar wind) 


16000 


510 


240 


490 



The velocity distribution of H atoms fn{r,W}i,t) may be calculated from the linear 
kinetic equation introduced in [ [5T|j : 

<9/h . _ df H F dfu f,^ .|HP f /^ u - m 

+f P {r, wh) J |wh - w u \afJf H (f, w^)dw^ - (u ph + u impiiCt ) f H (f, w H ). 

Here /h^uJh) is the distribution function of H atoms; f p (r,w p ) is the local distribution 
function of protons; w p and wu are the individual proton and H atom velocities, respec- 
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tively; cr^, is the charge exchange cross section of an H atom with a proton; v p h is the 
photoionization rate; m H is the atomic mass; z/ impact is the electron impact ionization 



force. The plasma and neutral components interact mainly by charge exchange. However, 
photoionization, solar gravitation, and radiation pressure, which are taken into account 
in equation ([!]), are important at small heliocentric distances. Electron impact ionization 
may be important in the heliosheath (region 2). The interaction of the plasma and H 
atom components leads to the mutual exchanges of mass, momentum and energy. These 
exchanges should be taken into account in the plasma equations through source terms, 
which are integrals of fn {r, w, t) . 

3.2. Solar wind and interstellar electron and proton components 

Basic assumptions necessary to employ a hydrodynamic approach for space plasmas 
were reviewed in |nj. In particular, it was concluded in the paper that interstellar and 
solar wind plasmas can be treated hydrodynamically. Indeed, the mean free path of the 
charged particles in the local interstellar plasma is less than 1 AU, which is much smaller 
than the size of the heliospheric interface itself. Therefore, the local interstellar plasma is 
collisional plasma, and a hydrodynamic approach can be used to describe it. Solar wind 
plasma is collisionless, because the mean free path of the solar wind particles is much 
larger than the size of the heliopause. Therefore, the heliospheric termination shock (TS) 
is a collisionless shock. A hydrodynamic approach can be justified for collisionless plasmas 
when scattering of charged particles on plasma fluctuations is efficient (" collective plasma 
processes"). In this case, the mean free path I with respect to collisions is replaced by l co a, 
the mean free path of collective processes, which is assumed to be less than the character- 
istic length of the problem L: l co u << L. However, the integral of "collective collisions" is 
too complicated to be used to calculate the transport coefficient for collisionless plasmas. 

One-fluid description of heliospheric and interstellar plasmas is commonly used in the 
global models of the heliospheric interface. However, since measurements of the solar wind 
show different electron and proton temperatures, two-fluid approach is more appropriate. 
The temperatures may remain different up to the termination shock and beyond due to 
the weak energy exchange between protons and electrons. 

Hydrodynamic Euler equations for proton and electron components, which take into 
account the influence of other components such as interstellar H atoms, pickup ions, 
cosmic rays, electric and magnetic fields, are written below. Mass balance or continuity 
equations are 



Index e denotes electrons, index p denotes solar wind protons. qi >e = uh ■ (v ph + b> impact ), 
qi tP = — J ua^ (u) f p (w) fH(wH)dwdwH are sources and sinks due to charge exchange, pho- 
toionization and electron impact ionization. Here, u = \wh~ w | is the relative atom-proton 
velocity, and wh and w are individual velocities of H atoms and protons, respectively. Mo- 
mentum balance equations are 



rate; and F is the sum of the solar gravitational force and the solar radiation pressure 




(2) 




+ VP S + m s V • {n s V s ® V.) - n s e s (E + -[V s x B]) + £ R sr = m s g 2>s (3) 



dt 



r 
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where <f 2 , e = J(vph+Vimpact)w H fH(w H )dwH, <?2, P = - / 1 uafJ(u)w p f H (w H )fp(wp)dwHdWp 
and R sr is the rate of momentum transfer from the r-gas component to the s-gas com- 
ponent. The symbol ® represents the dyadic product. The momentum transfer term R sr 
can be expressed in a general form through the collision integral S sr of kinetic equation 
of the s-gas component: R sr = — J m s c s S sr dc s , where c s = w s — V s . That R sr + R rs = 
is a consequence of this definition of R sr . 

Heat balance equations have the following form: 

d f 3 \ f 3 \ 

di \2 Ps ) + V (.2^7 + P * V ' ^ = E Qsr + m s q 3>s - m s q %s ■ V s (4) 

-*2 

with g 3i6 = / (v ph + v impact ) ^f H {w H )dw H , 

-*2 

g 3)P = — J J Mcr^f (u)^y fuiwu) fp(wp)dw p dvH ■ Q sr is the heat source due to interactions 
between between particles of s and r components. Those terms can be expressed in a 

2 

general form through the collision term S sr of the kinetic equation Q sr = — J ^^■S , pr (ic s . 
As a result of this interaction, we have a relation connecting Q sr and Q rs : Q sr + Q rs = 
-R sr -(V s -V r )(Vs^r). 

System (0)-(§) should be added by the state equations: P a = n a kT a (a = e,p), where 
k is Boltzman constant, and Maxwell equations: 

ldB - 4-k- 
Vx£ = — ;\7 ■ E = 4irp e ;V x B = — j; V • 5 = (5) 

c at c 

where p e is the charge density, e is the charge of electron, and j is current density. 
The displacement current has been dropped. Note, that charge and current densities 
of all charged populations should be taken into account in (5). Neglecting cosmic ray 
charges and currents we have p e = e{n p + n pui — n e ) and j = e(n p V p — n e V e + np Ui Vp U i). 
The number density of pickup ions, n pui , and the bulk velocity of pickup ions, V pui are 
integrals of pickup proton velocity distribution function: n pui = J f pui (w)dw, V pui = 
(J ' wfp Ui {w)dw)/n pui . 

Note that in equations (0)-(f§) we assume that pickup electrons are indistinguishable 
from original solar wind electrons, while pickup protons are considered as a separate 
population. 

The expressions for various interaction terms R sr , Q s (s = e,p)must be specified. 
Electron-proton collision terms can be taken in the form given by Braginski |]22|| : 

Rep = U e (6) 

3tz tyl 

Qpe Qep Rep^e k(T e Tpj (7) 

T e m p 

Here, n e is the electron number density; T e and T p are electron and proton densities, 
respectively; m e and m p are the electron and proton masses; u e is the electron velocity 
relative to the proton rest frame. Parameter r e characterizes the coupling between elec- 
trons and protons and corresponds to the electron collision time for collisional plasma 



22]. In collisionless heliospheric plasma, additional assumptions are needed to determine 



r e ; otherwise, it can be considered as free parameter. 
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3.3. Pickup ions 

To study pickup ion dynamical influence on the distant solar wind, the termination 
shock structure and, finally, on the global heliospheric interface structure, details of the 
process of charged particle assimilation into the magnetized plasma are needed. A newly 
created ion under the influence of the steady solar wind electric and magnetic fields 
executes a cycloidal trajectory with the guiding center, which is drifting at the bulk 
velocity of the solar wind. Assuming that the gyroradius is much smaller than the typical 
scale length, one can average velocity distribution function over the gyratory motion. 
Initial ring-beam distribution of pickup ions is unstable. Basic processes that determine 
evolution of pickup ion distribution are pitch-angle scattering, energy diffusion in the 
wave field generated by both pickup ions and the solar wind waves, convection, adiabatic 
cooling in the expanding solar wind, and injection of newly ionized particles. The most 
general form of the relevant transport equation to describe the evolution of gyrotropic 
velocity distribution function f pui = f pui (t,r,v, /x) of pickup ions in a background plasma 
moving at a velocity V sw were written in |p3[] , p4fl . f pui is a function of the modulus of 
velocity in the solar wind rest frame, and fi is the cosine of pitch angle. 

Complete assimilation of pickup ions into the solar wind would result in a great in- 
crease in the temperature with increasing heliocentric distance, which is not observed. 
Therefore, the solar wind and pickup protons represent two distinct proton populations. 
Nevertheless, the radial temperature profile of protons measured by Voyager 2 shows a 
smaller decrease as compared with the adiabatic cooling. A fraction of heating of so- 
lar wind protons may be connected with pickup generated waves |[25|| . Many aspects of 
pickup ion evolution were studied ( e.g., |2l| ; for review, see |TIJ, ||). However, today it 



still seems to be impossible to take into account all details of the assimilation process of 
pickup ions into the solar wind in the global models of the heliospheric interface structure. 
Instead, one may try to use the hydrodynamic approach. In this approach, equations (|2])- 
dU) written for pickup ions represent the balance of their mass, momentum and energy. 
The right sides of the equations include sources of pickup ions due to ionization processes: 

qi, P ui = n H u ph + J uafJ(u)f H (w H )f p (w)dwdv H 

q2, P ui = J {v P h + v impact )v H f H {v H )dv H + J J uafj (u)w H f H (w H ) f p (w p )dw H dw p + 

{Vph + "impact) ^fH{w H )fp{Wp)dw p dw H 
p p r \jj\ TV ^ 

+ J I ua e ^(u)^—^ — L fH(w H )fpui(wi)dwidwH 

To complete the model, one should also specify interaction terms R pu i, r , Q P ui,r (r ^ pui). 
The specification of these terms for pickup ion-proton interactions requires analysis of the 
pickup process in detail at the kinetic level. Global models usually assume immediate 
assimilation of pickup ions into the solar wind (one-fluid model) or perfect co-moving of 
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these populations V p = V pu i and no exchange of energy Q pu i, P = (two- or three-fluid 
models). Energy exchange term of pickup with ACRs, Qpui,acr — —Qacr, P ui, is specified in 
next subsection. 

3.4. Cosmic rays 

The cosmic rays are coupled to background flow via scattering with plasma waves. The 
net effect is that the cosmic rays tend to be convected along with the background plasma 
as they diffuse through the magnetic irregularities carried by the background plasma. 
Both galactic and anomalous cosmic rays can be treated as populations with negligible 
mass density and sufficient energy density. At a hydrodynamical level, the cosmic rays 
may modify the wind flow via their pressure gradient VP C with the net energy transfer 
rate from fluid to the cosmic rays given by V ■ VP C . P c (r, t) — ^ J °° f c (r,p, t)wp 3 dp is a 
cosmic ray pressure; f c (r,p,t) is the isotropic velocity distribution of cosmic rays. 

The transport equation of these particles has the following form ||]: 

Here p is the modulus of the momentum of the particle; D is the diffusion coefficient 
in momentum space, often assumed to be zero; k is the tensor of spatial diffusion; V = 
U + Vdrift is the convection velocity; U is the plasma bulk velocity; Vdrijt is a drift velocity 
in the heliospheric or interstellar magnetic field; and S(f,p,t) is the source term. 

At the hydrodynamic level, the transport equation of the cosmic rays in the heliospheric 
interface is: 

dP 

—1 = V[WP C - lc {U + U dr )P c ] + ( 7c - l)U ■ VP C + Q acr , pui (r, t) (9) 

Here we assume that D = 0; Ud r is momentum-averaged drift velocity; 7 is the polytropic 
index; and Q a cr, P ui is the energy injection rate describing energy gains of the ACRs from 



pickup ions. Chalov and Fahr ( |p6fl , [[27]]) suggested that Q a cr, P ui = —otp pu idivU ', where a 



is a constant injection efficiency defined by the specific plasma properties |p7[| . a is set to 
zero for GCRs since no injection occurs into the GCR component. 

4. OVERVIEW OF HELIOSPHERIC INTERFACE MODELS 

Together with Maxwellian equations @ the Boltzman equation (|XJ) for interstellar H 
atoms; sets of hydrodynamic equations (0)-(f§) written for solar protons, electrons and 
pickup ions; and equation (Q) written for anomalous and cosmic ray components form 
a closed system of equations, when interaction terms R sr , Q sr are specified. Possible 
specification of the interaction terms is given in equations @, (|7]). We have to note here, 
that such a complete model has yet to be developed. However, in recent years, several 
groups have focused their efforts on theory and modeling in order to understand some 
effects separately from others. In particular, the influence of the interstellar magnetic field 
on the interface structure was studied in p8| - [pT| for the two-dimensional case and in 



and [33f for the three-dimensional case. Both interstellar and interplanetary magnetic 



fields were considered in [33] and |35|. A comparison of these MHD models was given 



recently in |36]. Latitudinal variations of the solar wind have been considered in ]37||. The 
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influence of the solar cycle variations on the heliospheric interface was studied in the 2D 
[ 42| and in for the 3D case. In spite of many interesting findings in the 



case m 
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papers cited above, these theoretical studies did not take into account the interstellar H 
atoms, or took them into account but under greatly simplified assumptions, as it was done 
in [34 , where velocity and temperature of interstellar H atoms were assumed as constants 



in the entire interface. 

Since most of the observational information on the heliospheric interface is connected 
with interstellar neutrals and their derivatives as pickup ions and ACRs, we will focus on 
the models, which include interstellar neutrals in a more appropriate way. These models 
can be separated into two types. Models of the first type (Table 5) use a simplified fluid 
(or multi-fluid) approach for interstellar H atoms. A kinetic approach was used in the 
models of the second type. Development of the fluid (or multi-fluid) models of H atoms 
was connected with the fact that fluid (or multi-fluid) approach is simpler for numerical 
realization. At the same time such an approach can lead to nonphysical results. Results 
of one of the most sophisticated multi-fluid models [4"3" were compared with the kinetic 
Baranov-Malama model in [[44]]. The comparison shows qualitative and quantitative dis- 
agreements in distributions of H atoms. At the same time, it was concluded in [|J| that 
the two models agreed on the distances to the termination shock, heliopause and bow 
shock in upwind, but not in positions of the termination shock in downwind. 

4.1. One-fluid plasma models 

One of the common features in the models 
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H - El 



0-@ 

is that proton, electron and pickup ion components were considered as one fluid. The 
great advantage of this approach is that its equations are considerably simpler than the 
three- and two-fluid approaches (see next subsection). A key assumption of this approach 
is immediate assimilation of pickup protons into the original solar protons. In other 
words, it is assumed that immediately after ionization one cannot distinguish between 
original solar protons and pickup protons. Another important assumption is that electron 
and proton components have equal temperatures, T e = T p . For quasineutral plasma 



(n p + n pui = n e + o(n e )) this means that the pressure of the electrons is equal to half of 
total pressure (P = n e kT e + {n p + n pui )kT p « 2n e kT e = 2P e ). Let us denote total density 



m P n. 



+ m p (n p + n pui ), bulk velocity V = (£ s m s n s V s )/p (s = e,p } pui). 



Governing equations for the one-fluid approach can be obtained by summarizing equa- 
tions (0)-© for indexes s = e,p,pui. Introducing solar wind protons and pickup ions 



as co-moving and taking into account that m e 
general form: 



<< m p yield one-fluid equations in their 



Ql, Ql — fn p nH(^ph + ^impact) 



djpV) 

dt 



+ VP + V ■ (pV (g) V) - p e E [j x B) = q 2 - V ■ n e m e u e u e 



(10) 



Here j = £ Q n a e a V a , p e = J2 a n a e a ; q 2 = £ s m s q 2 ^ u s = V s - V. 

Wt il p ) + v ■ (l p9 ) ~ 9 ■ VP = " v (l PA ) + 93 " * ' 9 + 



(12) 



to 



Table 5: Models with multi-fluid approaches for interstellar H atoms. 



Reference 



GCR 



ACR IMF 



HMF Latitud. SW Time Pickup and 
asymmetry Depend. SW protons 



H atoms 



MULTI-FLUID 
Liewer et al., 1995 [46] 
Zank et al, 1996 [43] 
Pauls and Zank, 1997 [37] 
McNutt et al., 1998, 1999 [47], [4* 
Wang and Belcher, 1999 [42] 
Fahr et al, 2000 [49] 



APPROACH FOR H ATOMS 



+ 



+ 



+ 
+ 



+ 

* 
+ 



+ 



+ 



KINETIC 

Osterbart and Fahr, 1992 [50] 
Baranov and Malama, 1993 [51] 
Muller et al., 2000 [52] 



APPROACH FOR H ATOMS 



+ 



MyEshikov et al 



, 2^0^68] 
Aleksashov et CJ 2000 [69] 

□ n n 

Zaitsev anijHpmodenov, 2001 [70] 



+ 



+ 



+ 



one-fluid 
one-fluid 
one-fluid 
one-fluid 
one-fluid 
two- fluid 



one- fluid 
three-fluid 
one-fluid 
one- fluid 
one- fluid 
one- fluid 



No pickup ions not self- 
consistent 



one- fluid 



One-fluid 



one- fluid 



one-fluid 



one-fluid 



Monte Carlo 
with splitting 
particle 
mesh code 
Monte Carlo 
with splitting 
Monte Carlo 
with splitting 
Monte Carlo 
with splitting 
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> 1 — » — » ffi 

(j - Pe"^) • ( -E + -K X B + — (<?2,p + 

c en e 



where 53 = J2 S m sQ3,s- The relative electron velocity is connected with vector j as follows: 
j = p e V — en e u e . Note that to derive ([TJ) we use a generalized form of Ohm's law and 
neglect the terms proportional m e /m p in it: 



n e e (e + -V e x B^j = -VP e + R ep + m e (q 2)P + q% >pui - <f 2 , e 



If we assume momentum transfer term R ep as in (ffi) and neglect the term p e V in quasineu- 
tral plasma, Ohm's law may be rewritten: 

-* 1 — » — * Tfl — * 1 / 1 — * — » \ Tfl 

E = —[VxB] + e —j + [-j x B - VP e + — (g 2iP + g 2 , pui - g 2 , e ) (13) 

c r e n e e z en e \c J en e 

To derive the classical system of hydrodynamic equations applied for heliospheric in- 
terface in one-fluid models, one needs to ignore terms containing magnetic and electric 
fields in equations (PD and flTjp. If we also disregard the second-order term Vn e m e u e u e 
in flnp and V(5p e -u e /2) in (Tl2|), equations (|10|)-(|l2|) together with kinetic equation ([1|) for 



H atoms form a closed system of equations. 

To get from (|II])-(|IB|) a system of ideally conducting MHD equations with a magnetic 
field frozen in plasma, one obviously needs to make other assumptions in addition to that 
of high conductivity, when R m » 1. R m = AttuV L/c 2 is magnetic Reynolds number, 
with the electrical conductivity a = n e e 2 T e /m e , and V and L characteristic velocity and 
length, respectively. Vanishing electron pressure gradients, the Hall term, j x B and 
last term of ( |13D connected with the charge-exchange effect, corresponding terms in the 
generalized Ohm's law (|13|) also must be ignored. Under these conditions, the Ohm's law 
has its classical form E = —-[V x B], and last term of heat flux equation (|T2|) is equal 
to j 2 /a and can be neglected. Ideal MHD equations with source terms q±, 52,53, on the 
right-hand sides were considered in ||34j| , P7L |f48| , 



4.2. Three- and two-fluid plasma models 

For solar wind, the one-fluid model assumes essentially that wave-particle interactions 
are sufficient for pickup ions to assimilate quickly into the solar wind, becoming indistin- 
guishable from solar wind protons. However, as discussed above, Voyager observations 
have shown that this is probably not the case. Pickup ions are unlikely to be assimi- 
lated completely. Instead, two co-moving thermal populations can be expected. A model 
that distinguishes the pickup ions from the solar wind ions was suggested by Isenberg 
in ||71|| . Electrons were considered as a third fluid. The key assumption in the model is 



that pickup ions and solar wind protons are co- moving (V p = V pu i). It was also assumed 
that there is no exchange of thermal energy between solar wind protons and pickup ions. 
Isenberg's approach consists of two continuity equations (0) for solar protons and pickup 
ions; one momentum equation (|Tl|) and three energy equations (f|) for solar wind protons, 
electrons and pickup ions. In (|TT|) Isenberg neglect term p e E, the last term and assumes 
that j = cV x B/(4n). In energy equations he disregards the energy exchange terms Q sr . 



Note that Isenberg used the simplified form of source terms suggested in [72] and applied 



these equations to the spherically symmetric solar wind upstream the termination shock. 
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Another two-fluid approach to model solar wind protons and pickup protons was de- 
veloped recently in j49|. This model also assumes that convection speed of pickup ions 
is identical to that of solar wind protons. The pressure of pickup ions is calculated by 
assuming a rectangular shape of the pickup ion isotropic distribution function. In this 
case, the pressure can be expressed through pickup ion density pi and solar wind bulk 
velocity V sw as 

P P m = P P mV s 2 j5. (14) 

Therefore, governing plasma equations are i) one-fluid equations for the mixture of solar 
protons, pickup ions and electrons; ii) continuity equation for pickup ions; hi) two trans- 
port equations for ACRs and GCRs. The influence of cosmic ray components was taken 
as terms -V{P AC r + Pgcr) and -V ■ V(P A cr + Pgcr) ~ aP pui div(V) in the right-hand 
side of the momentum and energy equations, respectively. 

5. BARANOV-MALAMA MODEL OF THE HELIOSPHERIC INTERFACE 

The first self-consistent model of the two-component (plasma and H atoms) LIC inter- 
action with the solar wind was developed by Baranov and Malama ||51|| . The interstellar 
wind is assumed to have uniform parallel flow in the model. The solar wind is assumed to 
be spherically symmetric at the Earth's orbit. Under these assumptions, the heliospheric 
interface has axisymmetric structure. 

Plasma and neutral components interact mainly by charge exchange. However, pho- 
toionization, solar gravity and solar radiation pressure, which are especially important in 
the vicinity of the Sun, are also taken into account. 

Kinetic and hydrodynamic approaches were used for the neutral and plasma compo- 
nents, respectively. The kinetic equation (|I] for neutrals is solved together with the Euler 
equations for one-fluid plasma (|10D - (|12|). The influence of the interstellar neutrals is 
taken into account in the right-hand side of the Euler equations that contain source terms 
qi, g*2, <?3, which are integrals of the H atom distribution function /#(Vh) and can be cal- 



culated directly by the Monte Carlo method | 5q| . The set of kinetic and Euler equations 
is solved by iterative procedure, as suggested in |p5| . Supersonic boundary conditions 
were used for the unperturbed interstellar plasma and for the solar wind plasma at the 
Earth's orbit. The velocity distribution of interstellar atoms is assumed to be Maxwellian 
in the unperturbed LIC. The model results are discussed below in this section. 

5.1. Plasma 

Interstellar atoms strongly influence the heliospheric interface structure. In the presence 
of interstellar neutrals, the heliospheric interface is much closer to the Sun than in a pure 
gas dynamical case (Figure 2). The termination shock becomes more spherical. The Mach 
disk and the complicated shock structure in the tail disappear. 

The supersonic plasma flows upstream of the bow and termination shocks are disturbed. 
The supersonic solar wind is disturbed by charge exchange with the interstellar neutrals. 
The new ions created by charge exchange are picked up by the solar wind magnetic field. 
The Baranov-Malama model assumes immediate assimilation of pickup ions into the solar 
wind plasma. The solar wind protons and pickup ions are treated as one-fluid, called the 
solar wind. The number density, velocity, temperature, and Mach number of the solar 
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Z (AU) Z (AU) 



Figure 2. Effect of the interstellar neutrals on the size and structure of the interface struc- 
ture, (a) The heliospheric interface pattern in the case of fully ionized local interstellar 
cloud (LIC), (b) the case of partly ionized LIC. BS is the bow shock. HP is the heliopause. 
TS is the termination shock. MD is the Mach disk. TD is the tangential discontinuity 
and RS is the reflected shock. 



wind are shown in Figure 3 A. The effect of charge exchange on the solar wind is significant. 
By the time the solar wind flow reaches the termination shock, it is decelerated (15-30 
%), strongly heated (5-8 times) and mass loaded (20-50 %) by the pickup ion component. 

The interstellar plasma flow is disturbed upstream of the bow shock by charge exchange 
with the secondary atoms originating in the solar wind and compressed interstellar plasma. 
Charge exchange results in the heating (40-70 %) and deceleration (15-30 %) of the 
interstellar plasma before it reaches the bow shock. The Mach number decreases and for 
a certain set of interstellar parameters {uh,lic » n PjLIC ) the bow shock may disappear. 
Solid curves on Figure 3B correspond to the small ionization degree of LIC [n v j (n p +n H ) = 
1/6). The bow shock almost disappears. 

The interstellar neutrals also modify the plasma structure in the heliosheath. In a pure 
gas dynamic case (without neutrals) the density and temperature of the postshock plasma 
are nearly constant. However, the charge exchange process leads to a large increase of 
the plasma number density and a decrease of its temperature (Figure 3C). The electron 
impact ionization process may influence the heliosheath plasma flow by increasing the 
gradient of the plasma density from the termination shock to the heliopause p4]| . The 
effects of interstellar atom influence on the heliosheath plasma flow may be important, in 
particular, for the interpretations of kHz radio emission detected by Voyager ([p7[, []55|1 ) 
and possible future heliospheric imaging in energetic neutral atom (ENA) fluxes 0. 

5.2. Atoms 

Charge exchange significantly disturbs the interstellar atom flow. Atoms newly created 
by charge exchange have velocities of their ion partners in charge exchange collisions. 
Therefore, the velocity distribution of these new atoms depends on the local plasma 
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Figure 3. Plasma density, velocity, temperature and Mach number upstream of the ter- 
mination shock (A), upstream of the bow shock (B), and in the heliosheath (C). The 
distributions are shown for the upwind direction. Solid curves correspond to nH,Lic=^-^ 
cm" 3 , n Pt Lic=0.04: cm -3 . Dashed curves correspond to riff,L/c=0.14 cm -3 , n Pt Lic— 0.10 
cm" 3 . V LIC =25.6 km/s, T L/C =7000 K 
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Figure 4. Number densities and veloci- 
ties of 4 atom populations as functions of 
heliocentric distance in the upwind direc- 
tion. 1 designates atoms created in the su- 
personic solar wind, 2 atoms created in the 
heliosheath, 3 atoms created in the dis- 
turbed interstellar plasma, and 4 original 
(or primary) interstellar atoms. Number 
densities are normalized to nu,LiCi veloci- 
ties are normalized to Vuc- It is assumed 
that nu.Lic — 0.2 cm -3 , n P)LIC = 0.04 
cm -3 . 




Figure 5. Velocity distributions of 
four atom populations at the termination 
shock in the upwind direction. (A) pri- 
mary interstellar atoms, (B) secondary in- 
terstellar atoms, (C) atoms created in the 
heliosheath, (D) atoms created in the su- 
personic solar wind. V z is the projection 
of velocity on the axis parallel to the LIC 
velocity vector. Negative values of V z in- 
dicate approach to the Sun. V x is the 
radial component of the projection of ve- 
locity vector on the perpendicular plane. 
V z , V x are in cm/sec. It is assumed that 
n H ,Lic = 0.2 cm -3 , n pM c = 0.04 cm -3 . 
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properties. It is convenient to distinguish four different populations of atoms depending 
on where in the heliospheric interface they originated. Population 1 is the atoms created 
in the supersonic solar wind. Population 2 is the atoms originating in the heliosheath. 
Population 3 is the atoms created in the disturbed interstellar wind. We will call original 
(or primary) interstellar atoms population 4. The number densities and mean velocities 
of these populations are shown in Figure 4 as the function of the heliocentric distance. 
The velocity distribution function of interstellar atoms /jf(wfr, r) can be represented as 
a sum of the distribution functions of these populations: /# = fn,i + fn,2 + fn,3 + fu,i- 
The Monte Carlo method allows us to calculate these four distribution functions. The 
velocity distributions of the interstellar atoms in the 12 selected points in the heliospheric 
interface were presented in . For example, the velocity distributions at the termination 
shock in the upwind direction are shown in Figure 5. Note that velocity distributions of 
H atoms in the heliosphere were also presented in |52|| . However, different populations of 
H atoms cannot be considered separately in mesh particle simulations of H atoms |53| . 



Original (or primary) interstellar atoms are significantly filtered (i.e. their number 
density is reduced) before reaching the termination shock (Figure 4A). Since slow atoms 
have a smaller mean free path as compared with fast atoms, they undergo more charge 
exchange. This kinetic effect, called "selection", results in a deviation of the interstellar 
distribution function from Maxwellian ( Figure 5A). The selection also results in ~10 % 
increase of the primary atom mean velocity to the termination shock (Figure 4C). 

The secondary interstellar atoms are created in the disturbed interstellar medium 
by charge exchange of primary interstellar neutrals and protons decelerated by the bow 
shock. The secondary interstellar atoms collectively make up the "H wall", a density 
increase at the heliopause. The "H wall" has been predicted in |5S[ and detected toward 
a Cen 0. At the termination shock, the number density of the secondary neutrals is 
comparable to the number density of the primary interstellar atoms (Figure 4A, dashed 
curve). The relative abundances of the secondary and primary atoms entering the helio- 
sphere vary with degree of interstellar ionization. It has been shown in [17] that the rela- 
tive abundance of the secondary interstellar atoms inside the termination shock increases 
with increasing interstellar proton number density. The bulk velocity of the population 
3 is about -18 -19 km/s. The sign "-" means that the population approaches the Sun. 
One can see that the velocity distribution of this population is not Maxwellian (Figure 
5B). The reason for the abrupt behavior of the velocity distribution for V z > is that 
the particles with significant positive V z velocities can reach the termination shock only 
from the downwind direction. The velocity distributions of different populations of H 



atoms were calculated in |66| for different directions from upwind. The fine structures of 
the velocity distribution of the primary and secondary interstellar populations vary with 
direction. These variations of the velocity distributions reflect the geometrical pattern of 
the heliospheric interface. The velocity distributions of the interstellar atoms can be a 
good diagnostics of the global structure of the heliospheric interface. 

The third population of the heliospheric neutrals is the neutrals created in the 
heliosheath from hot and compressed solar wind protons. The number density of this 
population is an order of magnitude smaller than the number densities of the primary and 
secondary interstellar atoms. This population has a minor importance for interpretations 
of Ly a and pickup ion measurements inside the heliosphere. However, some of these 
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atoms may probably be detected by Ly a hydrogen cell experiments due to their large 
Doppler shifts. Due to their high energies, the particles influence the plasma distributions 
in the LIC. Inside the termination shock the atoms propagate freely. Thus, these atoms 
can be the source of information on the plasma properties in the place of their birth, i.e. 
the heliosheath ||. 

The last population of heliospheric atoms is the atoms created in the supersonic 
solar wind. The number density of this atom population has a maximum at ~5 AU. 
At this distance, the number density of population 1 is about two orders of magnitude 
smaller than the number density of the interstellar atoms. Outside the termination shock 
the density decreases faster than 1/r 2 where r is the heliocentric distance (curve 1, Figure 
4B). The mean velocity of population 1 is about 450 km/sec, which corresponds to the 
bulk velocity of the supersonic solar wind. The velocity distribution of this population 
is not Maxwellian either (Figure 5D). The extended "tail" in the distribution function is 
caused by the solar wind plasma deceleration upstream of the termination shock. The 
"supersonic" atom population results in the plasma heating and deceleration upstream of 
the bow shock. This leads to the decrease of the Mach number ahead of the bow shock. 



5.3. Recent developments in the Baranov-Malama model 

The Baranov-Malama model, the basic results of which were discussed above, takes into 
account essentially two interstellar components: H atoms and charged particles. To apply 
this model to space experiments, one needs to evaluate how other possible components 
of the interstellar medium influence the results of this two-component model. Recently, 
several effects were taken into account in the frame of this axisymmetric model. 

The influence of the galactic cosmic rays on the heliospheric interface structure was 
studied recently in ||67|| , ||68|| . The study was done in the frame of two-component (plasma 



and GCRs) and three-component (plasma, H atoms and GCRs) models. For the two- 
component case it was found that cosmic rays could considerably modify the shape and 
structure of the solar wind termination shock and the bow shock and change the positions 
of the heliopause and the bow shock. At the same time, for the three- component model 
it was shown that the GCR influence on the plasma flows is negligible as compared 
with the influence of H atoms. The exception is the bow shock, a structure that can 
be strongly modified by the cosmic rays. It was also found ( |)49|| ; Alexashov, private 
communication) that an anomalous component does not have a significant effect on the 
position of the termination shock. However, ACRs may significantly reduce compression 



at the termination shock |49[ 



Effects of the interstellar magnetic field on the plasma flow and on distribution of H 
atoms in the interface were studied in [^] in the case of magnetic field parallel to the 
relative Sun/LIC velocity vector. In this case, the model remains axisymmetric. It was 
shown that effects of the the interstellar magnetic field on the positions of the termination 
and bow shocks and the heliopause are significantly smaller as compared to model with 
no atoms P5fl . The calculations were performed with various Alfven Mach numbers in the 
undisturbed LIC. It was found that the bow shock straightens out with decreasing Alfven 
Mach number (increasing magnetic field strength in LIC). It approaches the Sun near the 
symmetry axis, but recedes from it on the flanks. By contrast, the nose of the heliopause 
recedes from the Sun due to tension of magnetic field lines, while the heliopause in its 
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wings approaches the Sun under magnetic pressure. As a result, the region of compressed 
interstellar medium around the heliopause (or "pileup region") decreases by almost 30 
%, as the magnetic field increases from zero to 3.5 xlO -6 Gauss. It was also shown 
in |)9|] that H atom filtration and heliospheric distributions of primary and secondary 
interstellar atoms are virtually unchanged over the entire assumed range of the interstellar 
magnetic field (0 - 3.5 ■ 10~ 6 Gauss). The magnetic field has the strongest effect on density 
distribution of population 2 of H atoms, which increases by a factor of almost 1.5 as the 
interstellar magnetic field increases from zero to 3.5 • 10~ 6 Gauss. 

Very recently a new non-stationary model of the solar wind interaction with two- 
component (H atoms and plasma) LIC was proposed in j70fl . In this model the primary 
and secondary interstellar atoms (populations 3 and 4) were treated as quasi-stationary 
kinetic gases. Population 1 of atoms originating in the supersonic solar wind was con- 
sidered as zero-pressure fluid. The calculations show that the qualitative features of the 
non-stationary SW/LIC interaction established in j|l| remain, but the effect of the solar 
activity cycle is quantitatively stronger because the interface is closer to the Sun than in 
the model with no atoms. The motion of the termination shock during the solar cycle on 
the axis of symmetry is about 30 AU. Due to the solar cycle variations of the neutralized 
solar wind (i.e. atoms of population 1) the region between the heliopause and the bow 
shock widens and the mean plasma density in the region becomes smaller than for the 
stationary problem. 



6. INTERPRETATIONS OF SPACECRAFT EXPERIMENTS ON THE BA- 
SIS OF THE BARANOV-MALAMA MODEL 

The Sun/LIC relative velocity and the LIC temperature are now well constrained ( [|i~3|l , 
Using the SWICS pickup He results and an interstellar HI/Hel ratio of 13 
± 1 (the average value of the ratio toward the nearby white dwarfs), Gloeckler et al.||19| 
concluded that n LIC (HI) = 0.2 ±0.03 cm ~ 3 . This estimation of Hue (HI) is independent 
of the heliospheric interface model but mo del- dependent for determination of the number 
density of H atoms from pickup fluxes. Estimates of interstellar electron number density 
require a theoretical model of the heliospheric interface. The Baranov-Malama model was 
used in |T7[] to study the sensitivity of the various types of indirect diagnostics of local 



interstellar plasma density. The diagnostics are the degree of filtration, the temperature 
and the velocity of the interstellar H atoms in the outer heliosphere (at the termination 
shock), the distances to the termination shock, the heliopause, and the bow shock, and 
the plasma frequencies in the LIC, at the bow shock and in the maximum compression 
region around the heliopause, which constitutes the "barrier" for radio waves formed 



in the interstellar medium. We also searched [17] for a number density of interstellar 
protons compatible with SWICS/Ulysses pickup ion observations, backscattered solar Ly 
a observed by SOHO, Voyager and HST, and kHz radiations observed by Voyager. Table 
1 presents the ranges of n Pi Lic obtained on the basis of the Baranov-Malama model and 
comparable to these observations. 



From analysis of the ranges, it was concluded in |17| that it is difficult in the frame of the 
model to reconcile the results obtained from all types of data as they stand now. There 
is a need for some modifications of the interpretations or of the confidence intervals. 
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Table 6 

Intervals of Possible Interstellar Proton Number Densities 



Type of Heliospheric Interface DiagnosticsRange of Interstellar Proton Number Density 



SWICS/Ulysses pick-up ion p| 
0.09 cm" 3 < n H ,TS < 0.14 cm" 3 
Ly- a, intensity [pf 



0.11 cm" 3 < u h ,ts < 0.17 cm" 3 
Ly-a, Doppler shift p§| - [|6l 



18 km s" 1 < V h ,ts < 21 km s" 1 
Voyager kHz emission (events) 

110 AU < i?Au < 160 AU 
Voyager kHz emission (cutoff) 

1.8 kHz 



0.02 cm" 3 < n Pi Lic < 0.1 cm" 3 
%>,lic < 0.04 cm" 3 or 
0.07 cm" 3 < 72 Pi lic < 0.2 cm" 3 
0.08 cm" 3 < n pMC < 0.22 cm" 3 
™p,lic = 0.04 cm" 3 



Two mutually exclusive solutions have been suggested: (1) It is possible to reconcile 
the pickup ions and Ly a measurements with the radio emission time delays if a small 
additional interstellar (magnetic or low-energy cosmic ray) pressure is added to the main 
plasma pressure. In this case, n Pi Lic = 0.07 cm" 3 and n Hi Lic — 0.23 cm" 3 is the favored 
pair of interstellar densities. However, in this case, the low frequency cutoff at 1.8 kHz 
does not correspond to the interstellar plasma density, and one has to search for another 
explanation. (2) The low-frequency cutoff at 1.8 kHz constrains the interstellar plasma 
density, i.e., n Pj Lic = 0.04 cm" 3 . In this case, the bulk velocity deduced from the Ly a 
spectral measurement is underestimated by about 30-50% (the deceleration is about 3 
km s" 1 instead of 5-6 km s" 1 ). Model limitations (e.g. a stationary hot model to derive 
the bulk velocity) or the influence of a strong solar Ly a radiation pressure may play a 
role. In this gnificant additional interstellar (magnetic or cosmic ray) pressure 

as compared with case (1) would be needed. 



This need for an additional pressure is in agreement with the conclusions made in [52 
which were derived from the analysis of the H wall absorption toward alpha Centauri . 
In their model the authors modified the equation of the state of the gas to simulate the 
effect of the interstellar magnetic field (IMF) and concluded that H wall absorption favors 
the "subsonic case". However, the best model of these authors corresponds to a neutral 
H density of 0.025 cm" 3 in the inner heliosphere, at least 4 times smaller than the density 
derived from the pickup ions. Also, the precision required to model the differences between 
the theoretical absorptions, namely small differences of the order of a few kilometers per 
second at the bottom of the lines, is of the order of the differences between the kinetic and 
multi-fluid model results for the same parameters in the supersonic case (see Appendix B 
in E5j ; f44|). Thus, an additional study of the absorption toward nearby stars for more 



realistic densities and models is desired. 



7. PROBLEMS FOR FUTURE WORK 



The Local Interstellar Medium interacts with the solar wind and influences the outer 
heliosphere in a complicated way. Several particle populations and magnetic fields are 
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involved in this interaction. From the interstellar side, the interacting populations are the 
plasma (electron and proton) component, H atom component, interstellar magnetic field, 
and galactic cosmic rays. Heliospheric plasma consists of original solar wind protons, 
electrons, pickup protons, and the anomalous component of cosmic rays. A large effort 
has been done to study the theoretical physics of the interaction region. However, a 
complete, self-consistent model of the heliospheric interface has not yet been constructed, 
because of the difficulty connecting both the multi-fluid nature of the heliosphere and 
the requirements of the different theoretical approaches for different components of the 
interaction. Many aspects were studied and reported here in previous sections. However, 
some aspects require additional theoretical explorations. Most theoretical models employ 
the one-fluid approach for solar wind and interstellar plasmas. It has been shown that, to 
derive one-fluid approach equations, several assumptions are needed. A key assumption 
that looks reasonable is co-moving character of all components. Another assumption for 
a one-fluid plasma model is the immediate assimilation of the pickup ion component into 
the solar wind. As demonstrated by space experiments, this is not the case and it would 
be more natural to consider solar protons and pickup protons separately as co-moving 
populations. The electron component should also be treated as a distinct population. 
However, since the assumption of the co-moving character of these three heliospheric 
plasma populations looks reasonable, the one-fluid approach gives us a reasonably accurate 
picture of the flow pattern (positions of the shocks and heliopause) and plasma velocity 
distributions. Theoretical models of pickup ion acceleration and diffusion can be employed 
to determine the distribution of thermal energy between solar wind and pickup proton 
components. A similar study should be done for electrons. 

Another important aspect of the solar/wind interaction is a study of the tail region 
of the solar wind and interstellar medium interaction. Although some studies were done 
( J76| , J77| ) it is still not clear at which heliocentric distances the gas (plasma and H atoms) 
parameters become indistinguishable from local interstellar parameters, or in other words, 
how far signatures of the solar system are noticeable in the interstellar medium. It is still 
not clear which of the two competiting processes is the most important in the tail region 
- charge exchange or plasma transport across the heliopause due to different instabilities. 
Studies of Saturn's and Earth's magnetic tails show that such tails can be very extended 

Finally, growing interest in heliospheric interface studies is connected with expectations 
that Voyager 1 will cross the termination shock soon. Many predictions of the time of the 
termination shock crossing by Voyager appeared in the literature. However, it seems that 
much more work should be done to explain and reconcile all available indirect observations 
of the heliospheric interface based on the unique model of the heliospheric interface. This 
work should be done especially because NASA plans to send a spacecraft to a heliocentric 
distance of at least 200 AU with a flight-time of only 10 or 15 years. Intensive theoretical 
study will help to optimize goals, instrumentation, and, finally, the scientific profit of this 
" interstellar" mission. 
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